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ABSTRACT 

The  rigid  body  motion  of  an  ellipsoid  in  a  second  order  fluid  (SOF)  under 
the  action  of  specified  (time  independent)  external  forces  and  torques  have 
been  obtained  to  first  order  in  the  Weissenberg  number  by  inverting  the 
resistance  relations  for  the  force  and  torque  under  specified  rigid  body 
motions.  The  reciprocal  theorem  of  Lorentz  was  used  to  bypass  the  calculation 
of  the  o(W)  velocity  field.  The  results  agree  with  known  analytic  solutions 
for  a  SOF  with  the  secondary  to  primary  normal  stress  ratio  of  “1/2. 

The  solution  procedure  was  also  tested  by  computing  the  torque  on  a 
translating  prolate  spheroid  with  aspect  ratios  ranging  from  slender  bodies  to 
near-spheres.  One  result  is  that  for  a  SOF  with  zero  secondary  normal  stress 
(Weissenberg  fluid),  previous  asymptotic  results  for  near-spheres  were  found 

to  be  accurate  even  at  fairly  large  aspect  ratios  (e.g.  2). 

New  results  for  non-degenerate  ellipsoids  suggest  that  the  orientation 
(as  monitored  by  Euler  angles)  and  trajectory  of  sedimenting,  non-axisymmetric 
particles  such  as  ellipsoids  provide  useful  information  on  the  rheology  of  the 
suspending  fluid.  / 

AMS  (MOS)  Subject  Classifications:  76D05,  35Q10 

Key  Words:  sedimentation,  spheroids,  ellipsoids,  viscoelastic  fluid, 
second-order  fluid 

Work  Unit  Number  2  -  Physical  Mathematics 


SIGNIFICANCE  AND  EXPLANATION 


Viscoelastic  fluids  are  encountered  in  many  manufacturing  processes , 

especially  in  the  form  of  solutions  and  melts  of  large  macro-molecules 
(polymers).  The  characterization  of  the  flow  properties  of  these  materials  is 

an  important  step  in  the  elucidation  of  the  molecular  structure  and 
processibility  of  the  material.  The  work  presented  here  show  how  data  from 

sedimentation  experiments  with  non-axisymmetric  particles  e.g.  ellipsoids,  may 
be  applied  to  determine  such  rheological  properties  (material  flow  properties) 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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1 .  Introduction 

It  has  been  shown,  both  theoretically  and  experimentally,  that  particles 
suspended  in  a  viscoelastic  fluid  behave  differently  than  when  suspended  in  a 
Newtonian  fluid  [1,2,3].  For  example,  an  orthotropic  particle  translating  in 
an  unbounded  Newtonian  fluid  does  not  experience  a  hydrodynamic  torque. 
Consequently,  such  particles  will  settle  in  a  Newtonian  fluid  without 
rotation.  On  the  other  hand,  experiments  have  confirmed  the  rotation  of 
needle-like  particle  settling  in  viscoelastic  fluids,  e.g.  Separan  solutions 
[1].  These  observations  are  consistent  with  theoretical  predictions  of  a 
torque  on  a  slender  particle  translating  in  a  second  order  fluid  (SOF).  This 
has  lead  earlier  investigators  to  consider  the  feasibilty  of  these  experiments 
as  rheometric  tools  for  dilute  polymeric  solutions. 

Here,  we  examine  the  motion  of  an  ellipsoidal  particle  in  a  viscoelastic 
fluid  in  the  limit  of  small  Welssenberg  number,  W.  The  problem  is  solved  by  a 
perturbation  procedure  in  which  the  zero-th  order  solution  is  the  Newtonian 
result  and  the  first  order  solution  includes  the  0(W)  effect  in  a  SOF.  The 
method  employed,  classified  as  the  "reciprocal  theorem  method"  by  Brunn  [3], 
bypasses  the  solution  of  the  0(W)  velocity  problem  and  requires  knowledge  of 
just  the  Newtonian  velocity  field. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-^-0041. 
material  is  based  upon  work  supported  by  the  National  Science  Foundation 
through  an  Engineering  Research  Initiation  Grant  CBT-8404451  and  the 
Presidential  Young  Investigator  Award  CBT-8451056,  and  DMS-8210950,  mod. 


The  solution  of  this  problem  is  an  extension  of  earlier  works  on  shapes 
such  as  spheres  and  needles.  However,  the  primary  objective  of  this  work  is 
to  determine  whether  motions  of  nonaxisymmetric  particles  can  differ 
sufficiently  from  that  of  axisymmetric  particles  to  the  extent  that  such 
measurements  could  be  usefully  applied  in  rheological  studies  of  the 
suspending  fluid.  To  this  end,  the  reader  must  keep  in  mind  that  the 
suspending  fluid  is  modelled  as  an  Incompressible  fluid  with  a  pre-specif led 
constitutive  equation  and  that  despite  similarities  in  the  techniques,  the 
present  work  is  not  an  effort  on  the  modelling  of  the  rheological  behavior  of 
suspensions. 

The  outline  of  the  paper  is  as  follows.  The  mathematical  formulation  of 
the  problem  is  presented  in  Section  2.  The  Newtonian  case  is  reviewed, 
including  recently  discovered  [4]  "singularity  solutions"  for  the  velocity 
field,  in  Section  3.  Results  which  have  been  obtained  for  the  non-Newtonian 
effects  are  described  in  Section  4,  including  the  effect  of  non-ax isymme try. 


i 


2.  Problem  Formulation 

In  the  limit  of  Theologically  alow  flows,  the  constitutive  equation  for  the 
suspending  fluid  reduces  to  the  n-th  order  Rivlin-Erickson  fluid  [5].  We 
shall  confine  our  attention  to  the  first  deviations  from  Newtonian  behavior  by 
approximating  the  fluid  as  a  second  order  fluid.  The  ultimate  test  of  whether 
such  a  regime  exists  in  real  polymeric  solutions  can  only  be  determined  by 
comparison  with  experiments.  The  bulk  stress  is  given  by  the  SOF  constitutive 
equation  [5] 

o  ■  -p6  +  pA,  +  (b. , -2b0)A. *A,  +  b0A0, 

where  p  is  the  viscosity,  b2  and  b1 1  are  material  properties  of  the  SOF  and 

the  A  are  the  Rivlin-Erickson  tensors  defined  by 
«*n 

^  -  Vv  ♦  (Vy)1,  A2  -  D^/Dt  «•  A^Vy)*  +  (Vv)‘Ar 

The  following  dimensional  analysis  defines  the  Weissenberg  number  for 
this  problem.  If  the  stresses  are  divided  by  the  viscous  scale,  pU/i,  where  U 
is  the  particle  translational  speed  and  l  is  a  characteristic  particle 
dimension,  then  the  constitutive  equation  becomes 

o  *  -p6  +  +  WA  -  t  +  WA, 

where  A  -  (1+a)A^*A^  -  ^A2 , 

and  W  -  (-2b2)U/(p£)  is  the  Weissenberg  number.  The  parameter  a  -  b11/(-2b2). 
The  definitions  for  W  and  a  are  motivated  by  the  fact  that  in  Couette  flow, 
the  SOF  has  the  primary  normal  stress  coefficient,  -  -2b2  and  the  secondary 
normal  stress  coefficient,  ¥2  ■  b^ ?  [5],  Experimental  results  with  polymeric 
systems  and  polymer  kinetic  theories  suggest  that  0  S  -a  S  0.2  C 3 » 5] . 


From  hereon,  the  A  are  dimensionless  and  defined  by 
«n 


Vu  +  (Vu) 


DA, /Dt  +  A, 

-I  “i 


(Vu)t  ♦ 


(Vu)*A 


1  • 


Throughout  this  paper,  we  will  use  v  to  denote  the  (dimensional)  velocity 
field  and  reserve  u  for  the  velocity  scaled  by  U. 

We  assume  that  inertial  effects  are  small,  and  also  assume  that  the  fluid 
is  weakly  elastic.  More  specifically,  we  require  Re  <<  W  <<  1 .  The  field 
variables  can  then  be  expanded  in  the  small  parameter  W: 


u-  u(0)  ♦  Wu(1)  ♦ 


.(0) 


Wp 


(1) 


and  substituted  into  the  equations  of  change  (momentum  and  mass  conservation), 
resulting  in  the  following  hierarchy  of  perturbation  problems: 


-vP(0)  +  v2u(0)  -  0,  V*u(0)  -  0, 


-Vp(1)  +  V2u(1)  -  -V-A|  (0),  V«u(1)  -  0. 


Thus  at  zero-th  order,  we  have  the  Newtonian  problem,  which  we  shall 
review  in  the  following  section,  while  the  first-order  problem  reduces  to  the 
forced  Stokes  equation  (which  will  be  considered  in  Section  4). 
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3.  The  Newtonian  Solutions 


Consider  an  ellipsoid  with  semlaxes  of  lengths  a,  b  and  c,  with  a  2  b  &  c  and 


We  choose  i,  -  a.  The  governing  equations  for  the  velocity,  y,  and  pressure, 
p,  are  the  Stokes  equations  and  the  equation  of  continuity  for  Incompressible 
flow, 


-Vp  +  yV^y  -  0, 

(3.2) 

V*  Y  -  0. 

(3-3) 

with  u  denoting  the  fluid  viscosity. 

Since  the  governing  equations  are  linear  in  this  case,  we  may  consider 
the  translational  and  rotational  motions  separately.  Furthermore,  the 
translational  and  rotational  problems  are  completely  decoupled  in  the  sense 
that  the  fluid  exerts  no  force  on  a  rotating  ellipsoid  and  no  torque  on  a 
translating  ellipsoid.  The  boundary  conditions  for  the  two  problems  are,  in 
the  Cartesian  coordinate  system  fixed  with  respect  to  the  particle  as 
specified  by  equation  (3.1): 

Translating  Ellipsoid 

v  -  0  on  the  particle  surface  (3.4a,b) 

y  -->  -U  as  | x|  — >  ®. 

Rotating  Ellipsoid 

y  -  0  on  the  surface  (3.5a,b) 

*  "_>  ■»"*  as  | s |  -->  ®. 
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The  force-translation  and  torque-rotation  relations  may  be  expressed  as 

[6,7,8], 

F  *  -ytK-U  and  T  -  -pr R*», 

P 

where  the  components  of  the  resistance  tensors  K  and  R  with  respect  to  the 

t  r* 

particle-fixed  axes  are  given  by  -  R^  -  0  for  i*j  and 


2  -1 

l6Trabc(x0  ♦  aQa  )  f 

(3.6) 

l|irabc(b280  +  c2Y0)_1 (b2*c2) . 

(3.7) 

Here,  xQ.  and  YQ  are  constants  which  are  obtained  by  evaluating  the 

following  ellipsoidal  harmonics  at  A-0.1 

X(A)  -  abc  [P(t)]_1dt,  (3-8) 

J  A 
f0» 

a(X)  -  abc  [(a2+t)P(t)]_1dt,  (3-9) 

'  A 

with  P(t)  -  [(a2+t)(b2+t)(c2+t)]1/2. 

The  lower  limit  of  the  definite  integral,  A(x,y,z),  is  the  ellipsoidal 
coordinate  which  is  the  positive  root  of 

x2  y2  z2 
a*+A  +  b2+A  *  c*+A  *  1  * 

The  functions  8(A),  Y(A)  and  the  expressions  for  the  other  diagonal  elements 

t  p 

of  K  and  R  are  obtained  by  successive  cycling  of  the  dependence  on  a,  b  and 
c  and  the  subscripts  1 ,  2  and  3- 


1.  a,  8  and  Y  are  as  defined  in  Happel  &  Brenner  [9]  and  differ  from 
Jeffery's  [8]  definition  by  a  factor  of  (abc). 


We  shall  see  in  the  next  section  that  the  Newtonian  velocity  field  is 
required  in  the  calculations  for  the  first  non-Newtonian  terms  in  the 
expansion  solution.  The  classical  solution  of  Oberbeck  [6]  for  a  translating 
ellipsoid  and  Edwardes  [7]  for  a  rotating  ellipsoid  are  rewritten  as  velocily 
fields  forced  by  a  distribution  of  Stokes  singularities  (hence  the  name 
"singularity  solutions")  in  order  to  highlight  the  common  structure  of  the  two 
solutions.  The  velocity  field  of  the  point  force,  i.e.,  the  fundamental 
solution  of  the  Stokes  equation,  is  I/(8iru),  where  I,  the  Oseen-Burgers  tensor 
defined  by 


[u  ‘  Kj  *  PVy  "lth 


X  • 


and  its  pressure  field. 


2y 

PJ  ‘  7"S  • 


satisfy  the  Stokes  equation  with  point  forcing, 


-  H{  *  “,2lu 

and  the  continuity  equation. 


8iru6^  j6(x) , 


■  °- 


(3.10a) 

(3.10b) 


(3.11a) 

(3.11b) 


as  shown  in  reference  [9],  Chapter  2. 

The  Newtonian  solutions  of  Oberbeck  and  Edwardes  can  then  be  rewritten  as 

C4D: 


v(x)  -  L(n)*||  f(n)(x',y')il  +  ^§^2}  I(x-x*)/( 8ttm)  dx’dy’,  (3-12) 

E 

with  n-1  for  translation  and  n-2  for  rotation.  The  two-dimensional 
Integration  is  performed  over  the  fundamental  ellipse,  E,  which  will  be 
defined  below.  The  terms  in  equation  (3.12)  are: 


i 

1 


l 


J 
J 
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"V 


and 


f  (XV),  (2n-_1_)  2n-3 

f(n)U,y;  2TraEbE  q 

Q(x.y)  -  [1  -£-£]1/2. 

aE  °E 


,  2  2,1/2 

aE  -  (a  -c  ) 


.  ..2  2,1/2 

t>E  -  (b  -c  )  , 


(n) 


C  ~  F  if  n-1 , 

1  S*V  ♦  T*7  if  n-2 


(3.13a-e) 


For  n»1 ,  L^*I  -  -F«I,  so  the  translational  solution  is  generated  by  a 

distribution  of  Stokes  monopoles  (or  Stokeslets)  over  the  fundamental  ellipse 

(2) 

with  density  f^j.  For  n=2,  L  *1  -  [(S+T)*7]*I,  so  the  rotational  solution 
is  generated  by  the  Stokes  dipole  field,  S+T  where  S  and  T  are  respectively, 
the  symmetric  (stresslet)  and  anti-symmetric  parts  of  the  stress  dipole, 

|  (o*n)x  dA. 

The  dyadic  T  is  related  to  T,  the  hydrodynamic  torque  exerted  on  the 
ellipsoid,  by 

Tij  “  “  2®iJkTk* 

The  domain  of  integration,  E(x',y'),  is  the  interior  of  the  fundamental 
ellipse, 


The  fundamental  ellipse  is  the  degenerate  elliptical  disk  in  the  family  of 
ellipsoids  that  are  confocal  to  the  particle  ellipsoid.  The  major  and  minor 
semi-axes  of  the  fundamental  ellipse,  aE  and  b£,  are  given  in  equation 
(3- 13c)  .  The  density  function  f(n  in  E(x',y’)  may  also  be 

2.  References  [10,11]  use  k  and  (k*-h2)1/2  in  lace  f  our 

E  ana  d_ 
c*  • 


1 


a 
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interpreted  as  the  surface  singularity  distribution  for  an  elliptical  disk,  as 

can  be  seen  by  looking  at  the  limit  c  -  0  in  equation  (3-12). 

The  function  q(x,y)  which  appears  in  f^  plays  an  analogous  role  in  the 

potential  theory  for  ellipsoidal  particles  (see  [10]).  In  fact,  in  potential 

theory,  q  1  is  the  requisite  charge  distribution  over  the  fundamental  ellipse 

which  generates  ellipsoidal  equipotential  surfaces.  Chwang  &  Wu  [12]  have 

noted  that  the  distribution  of  Stokes  multipoles  in  low-Reynolds-number 

problems  is  similar  to  the  distribution  of  multipoles  in  analogous  problems  in 

potential  theory,  except  for  the  presence  of  additional  degenerate  multipoles 
2 

(the  V  I  term)  in  equation  (3.12).  The  presence  of  such  quadrupoles  (or 
potential  doublet)  when  n-1  and  octupoles  when  n-2  in  (3-12)  are  consistent 
with  (and  in  fact  extend)  the  rules  stated  in  [12]  for  prolate  spheroids. 

To  complete  the  solution,  we  must  relate  F,  T  and  S  in  terms  of  the 
knowns,  0  and  u.  The  expressions  for  F  and  T  are  obtained  from  equations 
(3.6)  and  (3.7).  The  stresslet  is  identically  zero  for  a  sphere  undergoing 
rigid  body  motions  in  a  quiescent  fluid.  However,  for  a  nonspherical 
ellipsoid,  a  stresslet  is  generated  by  a  rotational  motion  [13].  The  linear 
relation  between  the  stresslet  and  the  angular  velocity  can  be  deduced  from 
[8]  as 

s  *  -urg*«, 

where  rQijk  *  0  unless  {i,j,k}  is  a  permutation  of  (1,2,3)  and 

rQ123  =  rQ213  -  -  |iruabc(a2a0  +  b280)_1  (a2-b2) .  (3.1*0 

The  remaining  elements  of  rg  are  obtained  by  cycling  1,2,3  and  a,b,c.  This 
completes  the  description  of  the  singularity  solution.  Further 
simplifications  of  these  shape-dependent  tensors  are  possible  in  the 
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degenerate  case  of  ellipsoids  of  revolution.  These  and  other  asymptotic 
expressions  for  ellipsoidal  and  spheroidal  resistance  functions  are  provided 
in  [U]. 


In  summary,  the  basic  results  for  the  Newtonian  velocity  field  are: 

1)  The  disturbance  velocity  field  for  a  translating  ellipsoid  (or  a  fixed 
ellipsoid  in  a  uniform  stream)  is  generated  by  a  distribution  of  3tokes- 
lets  and  potential  doublets  over  the  fundamental  ellipse. 

2)  The  disturbance  fields  for  a  rotating  ellipsoid  or  a  fixed  ellipsoid  in 
a  constant  vorticity  field  (and  for  a  stationary  ellipsoid  in  a  rate-of- 
strain  field)  are  generated  by  a  distribution  of  rotlets,  stresslets  and 
Stokes-octupoles  over  the  fundamental  ellipse. 

3)  For  prolate  spheroids,  the  fundamental  ellipse  degenerates  into  a  line 
segment  from  one  focal  point  to  the  other  and  the  singularity  solutions  of 
Chwang  &  Wu  [12]  are  recovered.  Furthermore,  if  the  spheroid  is  slender, 
the  Chwang  &  Wu  distributions  reduce  into  forms  that  are  analogous  to 
those  used  by  Leal  [1]. 


In  actual  computation,  e.g.  evaluation  of  the  volume  integrals  over  the  fluid 
region  in  the  following  section,  it  is  more  efficient  to  first  Integrate  these 
representations  over  the  fundamental  ellipse  to  obtain  the  classical  solutions 
(i.e.  use  the  classical  forms).  However,  in  the  present  form,  these 
representations  are  more  amenable  to  analyses  such  a3  asymptotic  expansions  in 
the  limit  of  degenerate  shapes.  The  present  forms  also  provide  a 
straightforward  bridge  to  the  earlier  works  on  slender  bodies. 


i 


i 
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4.  The  Non-Newtonian  Contributions 


4.1  The  Torque  on  a  Translating  Ellipsoid 

We  first  derive  the  expression  for  the  torque  on  a  translating  ellipsoid. 

This  relation  is  then  "inverted"  to  determine  the  rotational  motions  of  a 
torque-free,  sedimenting  ellipsoid. 

The  perturbation  expansion  for  the  torque  is: 

T  =*  +  WT*1*  ♦  ...  (4.1) 


As  discussed  in  the  previous  section,  the  Newtonian  term  is  identically  zero 
for  a  translating  ellipsoid.  The  expression  for  the  0(W)  term  is: 


uU*2{ 


JS  * 


xx(A*n)  dA  + 


n)  dA  }, 


(4.2) 


The  first  Integral  in  (4.2)  is  simply  the  contribution  from  the  non-Newtonian 
terms  in  the  constitutive  equation,  evaluated  with  the  Newtonian  (zero-th 
order)  velocity  field.  The  second  integral  is  the  contribution  from  the 
stress  field  of  Wu^.  Obviously,  we  need  to  retain  only  the  Newtonian  terms 
in  this  stress  field.  Thus,  it  would  appear  that,  given  the  Newtonian 
solution,  the  first  term  can  be  calculated  directly,  whereas  the  second  term 
would  require  the  solution  of  the  forced  Stokes  equation.  However,  as  shown 
in  [1],  the  Lorentz  reciprocal  theorem  may  be  applied  to  convert  the  second 
integral  into  an  expression  which  requires  only  the  Newtonian  velocity  field, 
thereby  bypassing  the  0(W)  velocity  problem.  Explicitly,  for  particles  of  any 
shape,  we  have: 


f 

Js 


eijkxj(5 


(D 


-)k 


dA 


J'r 


tikAki,t 


dV, 


(4.3a) 


where  t^  is  the  i-th  component  of  the  Newtonian  solution  for  the  particle 


rotating  with  unit  angular  velocity  in  the  j-th  coordinate  direction  and  the 
volume  integration  is  over  the  unbounded  fluid  region  exterior  to  the 
particle.  Brunn  [14]  has  derived  a  more  compact  formula  by  noting  that  if  tj^ 
is  used  in  equation  (4.3a)  instead  of  ,  an  integration  by  parts  cancels 
the  surface  integral  in  equation  (4.2).  Thus 


f  l(t 
JVf  2^ki,l 


t£i,k) 


AK£  dV* 


(4.3b) 


The  results  of  the  previous  section  imply  that  the  required  velocity  fields 
may  be  expressed  as: 


.(0) 


-u,  +  F, 


ff 

JJ 

E 


f(1)<x’ 


,y'){i 


^c2q2V2}  I1j(x-;')/(8tr)  dx’dy'. 


(4.4) 


fcik  ”  Q£jk+ 


lrR  ^ 
2  nkmemj£; 


(4.5) 


JJ  f(2)(x\y')fl^2q2V2}  I£t  J(;-x')/(8ir)dx*dy'. 

E 

The  resulting  torque-translation  relation  for  the  0(W)  term  may  now  be 
expressed  as 

T<in  ■  "IjkYk-  <"-6> 

as  can  be  seen  from  combining  equations  (4.2),  (4.3b)  and  (4.4).  Furthermore, 
n|jj  (no  sum)  is  identically  zero  for  i-1,2,3  because  of  particle  symmetry. 
Thus  the  0(W)  contribution  to  the  torque  on  a  translating  ellipsoid  in  a  SOF 
is  completely  characterized  by  three  shape-dependent  parameters,  23^N1 32^ ’ 

( ^231  ■*’N2i 3^  and  ^N312+N321)‘  furthermore,  the  N|jk  are  linear  with  respect  to 
a,  therefore,  the  complete  solution  may  be  determined  from  the  solution  at 
(only)  two  distinct  values  of  a.  Before  turning  to  the  sedimentation  problem, 
we  examine  the  special  case  when  a  -  -1/2,  for  which  an  analytical  solution  is 
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possible. 


Giesekus  [15]  has  noted  that  for  a  -  -1/2,  V»A  is  lrrotational  and  may  be 
written  as  the  gradient  of  a  scalar  potential.  The  0(W)  equation  may  be 
solved  directly  by  Incorporating  the  right-hand-side  into  the  pressure.  Brunn 
[16]  has  applied  this  simplification  to  the  general  motion  of  a  rigid  particle 
in  a  SOF.  For  a  particle  undergoing  steady  translation,  his  results  reduce  to 
the  extremely  simple  formula, 

T(1)  -  jU/U)  F(0)*y.  (4.7) 

For  ellipsoids,  this  reduces  to 

T<U  -  ^wa/UK^g-^H^Ug,  (4.8) 

and  two  analogous  expressions  obtained  by  cycling  the  Indices.  Necessarily, 
the  general  solution  described  below  must  reproduce  these  results  when 
o  ■  -1/2. 

Results  for  the  torque  on  a  translating  ellipsoid  were  obtained  by 
numerical  integration  of  the  right-hand-side  of  equation  (4.3b)  using  an 
adaptive  quadrature  package  provided  by  the  Numerical  Algorithms  Group.  These 
results  were  also  reproduced  using  three  nested  one-dimensional  Integrations 
in  the  ellipsoidal  coordinates,  using  the  intuitively  obvious  procedure  of 
finer  grid  sizes  in  the  region  near  the  surface.  The  linearity  with  respect 
to  a  was  exploited,  i.e.,  the  integrand  was  split  into  terms  independent  of 
and  first  power  in  o. 

The  (dimensionless)  torque  on  prolate  spheroids  with  shapes  ranging  from 
slender  bodies  to  near-spheres  is  presented  in  Figure  1 .  The  torque  for 
arbitrary  a  may  be  obtained  from  the  figure  by  using  the  plots  for  a  -  0  and 
the  a-coeff icient.  For  a  -  0,  the  asymptotic  solution  for  near-spheres 
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-V-V 


remains  a  good  approximation  even  at  b/a  -  0.5.  In  contrast,  there  is  a 
significant  deviation  of  the  a-coeff icient  from  the  asymptotic  solution  (the 
x-axis)  at  this  value  of  b/a.  When  o  »  -1/2,  the  results  agree  with  equation 
(4.8)  over  the  entire  range  of  a.  Finally,  the  -a  coefficient  vanishes  when 
b  -  0.  There  is  a  sharp  "boundary  layer"  in  the  curve  which  does  not  appear 
in  the  scale  of  Figure  1.  The  maximum  in  the  curve  lies  at  a  value  of  b/a 
between  0.01  and  0.001.  (In  comparison,  when  a  »  -1/2  the  maximum  is  near 
b/a  -0.15). 

Thus  the  results  presented  here  agree  with  those  in  reference  [14]  for 
near-spheres  and  satisfy  equation  (4.8)  over  the  entire  range  of  aspect 
ratios.  In  the  slender-body  limit,  the  results  do  not  approach  those 
presented  in  [1].  At  the  present  time,  I  suspect  that  this  disagreement  is 
due  to  an  error  in  [1],  in  part  because  that  result  does  not  agree  with  the 
analytic  solution  for  a  -  -1/2.  However,  it  should  be  noted  that  there  are 
differences  between  slender  spheroids  and  slender -bodies  generated  by  uniform 
distribution  of  stoke3lets  and  rotlets^. 

4.2  Sedimentation 

We  shall  first  consider  the  general  formulation  for  the  translational  and 
rotational  velocities  of  a  particle  whose  motion  is  driven  by  specified 
(time-independent)  forces  and  torques.  The  sedimentation  problem  is  a  special 
case.  The  relations  between  the  force  and  torque  and  the  translational  and 
rotational  velocities  are,  accurate  to  0(W), 

F1  ■  -“tKuuj  *  <HljkuA  *  mu*uA  *  *  0<w2)’  ("-9a) 

T.  Chwang  h  Wu  [17]  note  that  there  are  fundamental  differences  between 
slender  spheroids  and  those  generated  with  a  uniform  distribution  of 
rotlets. 


-u  Ry-J  -  (N1JkUjUk  ♦  N1JkUjUk  ♦  NjJkWjUk)W  ♦  0(W  ). 


(4.9b) 


T 


i 


For  ellipsoids  and  other  orthotropically-symmetric  particles,  M1 -0  as  shown  by 
[1].  The  consequences  for  sedimentation  are  now  considered. 

The  force-translation  relation  implies  an  0(1)  translational  velocity, 
which  when  inserted  into  equation  (4.9b)  implies  that  m  is  0(W)  smaller  than 
U.  Thus  (4.9a)  and  (4.9b)  may  be  Inverted  to  arrive  at  the  following 
expressions  for  the  sedimenation  and  rotational  velocities: 


Ui  *  (u^ij  )_1  FjXt  +  0(W*),  (4.10a) 

«1  -  (u^ij)’1  <NjkfcW  +  0(w2)-  (4.10b) 


Since  the  sedimentation  and  rotational  motions  are  small,  the  sedimentation 
trajectory  and  the  evolution  of  the  particle  orientatation  may  be  determined 
by  using  small  time  steps  and  (4.10a)  and  (4.10b).  The  trajectory  of  the 
ellipsoid  center  and  the  three  Euler  angles  (defined  by  Figure  2)  are  traced 
for  different  ratios  of  the  normal-stress  coefficients  (Figures  3a-3d)  for  the 
representative  case  with  b/a  -  0.25  and  c/a  -  0.05.  Several  points  pertaining 
o  the  use  of  non-axi symmetric  particles  are  raised  by  these  figures. 

The  projection  of  the  trajectory  of  the  particle  center  onto  a  horizontal 
plane  deserves  special  mention.  Since  the  ellipsoid  eventually  achieves  the 
"vertical"  orientation  and  falls  without  horizontal  drift,  the  plane  curve  of 
Figure  2  approaches  a  limit  point  at  large  t.  Moreover,  the  fluid  rheology 
affects  the  shape  of  this  plane  curve  only  through  the  ratio,  o.  In  contrast. 


for  axisymmetric  particles,  the  projected  curve  does  not  provide  as  much 


useful  Information  since  It  is  simply  a  straight  line  in  the  X-Y  plane.  A 
change  in  a  results  in  a  change  in  the  length,  but  not  the  direction  of  the 
line. 

Another  behavior  which  may  be  exploited  is  the  different  time  scales  in 
the  decay  of  the  Euler  angles  e  and  tp.  The  particle  of  Figure  3  represents 
the  case  where  41  decays  faster  than  6.  In  general,  the  two  time  scales  may  be 
adjusted  relative  to  each  other  by  varying  b/a  and  c/b. 

In  conclusion,  the  theoretical  analyses  presented  here  suggest  that  the 
sedimentation  configurations  and  trajectories  of  non-ax isymmetrlc  particles 
may  ion 

is  being  initiated  to  follow  this  lead. 
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Captions  for  Figures 

2  2 

1.  The  torque,  T  -  WpUa  (f  -  ag)U  U  /U  ,  on  a  prolate  spheroid 

z  x  y 

translating  in  a  SOF,  with  W  -  ^U/Cpa)  and  a  -  • 

f  ( — )  and  g  (-  -  -).  The  solid  curve  is  Brunn's  [1M]  asymptotic  result 
for  f(b/a)  for  near-spheres.  The  asymptotic  curve  for  g(b/a)  is  coincident 
with  the  horizontal  axis. 

2.  The  definition  for  the  Euler  angles  of  the  subsequent  figures.  The  body- 
fixed  axes,  x,  y,  and  z  are  obtained  from  the  laboratory-fixed  axes, 

X,  Y  and  Z  as  follows:  rotation  about  the  Z-axis  by  an  angle  4)  until 
the  X-axis  coincides  with  the  line-of -nodes,  followed  by  rotation  about 
the  line-of-nodes  by  an  angle  0  until  the  Z  axis  coincides  with  the 
x-axis.  The  third  Euler  angle,  4»t  is  the  angle  between  the  y-axis  and 
the  line-of-nodes. 

3.  The  initial  orientation  is  0-45° ,  <j>-0°  and  ip— ^45®  and  the 

initial  location  of  the  center  is  (0,0,0).  The  curves  are  for  a-0  ( - ), 

a--0.1  ( — )  and  a«-0.2  (-  -  -).  All  time  axes  have  been  scaled  with 
the  characteristic  time,  (a/U)W  1 . 

a.  The  trajectory  for  the  ellipsoid  center  projected  onto  the  X-Y  plane. 

b.  The  evolution  of  0. 

c.  The  evolution  of  #.  Without  loss  of  generality,  the  X  and  Y  axes  are 
defined  so  that  $-0  initially. 

d.  The  evolution  of  4>. 
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